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Abstract. - Systems described by n-component (f> 4 models in a cx> d_1 x L slab geometry of finite 
thickness L are considered at and above their bulk critical temperature T c ,oo- The renormalization-group 
improved perturbation theory commonly employed to investigate the fluctuation-induced forces ("thermo- 
dynamic Casimir effect") in d — 4 — e bulk dimensions is re-examined. It is found to be ill-defined beyond 
two-loop order because of infrared singularities when the boundary conditions are such that the free prop- 
agator in slab geometry involves a zero-energy mode at bulk criticality. This applies to periodic boundary 
conditions and the special-special ones corresponding to the critical enhancement of the surface interac- 
tions on both confining plates. The field theory is reorganized such that a small-e expansion results which 
remains well behaved down to T Ci00 . The leading contributions to the critical Casimir amplitudes A per 
and A sPjSp beyond two-loop order are ~ (u*)( 3-6 " 2 , where u* = 0(e) is the value of the renormalized 
4> 4 coupling at the infrared-stable fixed point. Besides integer powers of e, the small-e expansions of these 
amplitudes involve fractional powers e^ 2 , with k > 3, and powers of lne. Explicit results to order e 3 ^ 2 
are presented for A por and A sp , sp , which are used to estimate their values at d — 3. 



Fluctuations associated with long wave-length, low-energy excitations play a crucial role in deter- 
mining the physical properties of many macroscopic systems. When such fluctuations are confined by 
boundaries, walls, or size restrictions along one or several axes, important effective forces may result. 
In those cases where the continuous mode spectrum that emerges as the system becomes macroscopic 
in all directions is not separated from zero energy by a gap, these fluctuation-induced forces are long- 
ranged, decaying algebraically as a function of the relevant confinement length L (separation of walls, 
thickness of the system, etc). 

A prominent example of such forces are the Casimir forces [1] induced by vacuum fluctuations of 
the electromagnetic field between two metallic bodies a distance L apart [2-A]. Analogous long-range 
effective forces occur in condensed matter systems as the result of either (i) thermal fluctuations at 
continuous phase transitions or else (ii) Goldstone modes and similar "massless" excitations [5-12]. 
In particular the former ones, frequently called "critical Casimir forces," have attracted considerable 
theoretical and experimental attention recently. Beginning with the seminal paper by Fisher and de 
Gennes [5], they have been studied theoretically for more than a decade using renormalization group 
(RG) [6-10] and conformal field theory methods [13], exact solutions of models [12], as well as Monte 
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Carlo (MC) simulations [14-16]. Though their detailed experimental investigation just began, recent 
experiments [17-21] have provided clear evidence of their occurrence. 

The goal of the present work is to reexamine the existence of the e = 4 — d expansion for crit- 
ical Casimir forces. Considering n-component </> 4 models in oo d_1 x L slab geometries, subject to 
different kinds of large-scale boundary conditions (BC), we shall show that for some of these BC, 
the e expansion breaks down at the bulk critical temperature T CiOC beyond two-loop order because 
of infrared singularities. This breakdown occurs quite generally whenever the two-point correlation 
function in Landau theory has a zero mode at bulk criticality, as it does for periodic BC and for those 
corresponding to critical enhancements of the short-range surface interaction at both walls [22, 23]. 

As we shall show, the infrared problems one has in these cases require a special treatment of the 
zero mode. This leads to a modified RG-improved perturbation expansion, which is well-defined at 
and above T Cj00 , but yields contributions to the critical Casimir force of the form (u*)' 3 " 6 ^ 2 beyond 
two-loop order, where u* = 0(e) is the value of the renormalized </> 4 coupling constant at the infrared- 
stable fixed point in d = 4— e bulk dimensions. Thus, contrary to common belief, this quantity does not 
generally have an expansion in integer powers of e; the present theory requires fundamental revision. 

That zero modes may play an important role and require special treatment has been observed 
before in studies of finite-size effects in systems that are finite in all, or in all but one, directions [24]. 
Such systems (with short-range interactions) differ from the ones in slab geometry considered here in 
that they do not have a sharp phase transition at T > 0, except in the bulk limit. 

To put things in perspective, consider a system in a oo d_1 x L slab geometry whose bulk (L — > oo) 
critical behavior is representative of the universality class of the d-dimensional n-component </> 4 model 
with short-range interactions. Its reduced free energy per unit cross-sectional area A — > oo can be 
decomposed as 

^ = i im JiTt = L + f>* + f*> b + L ) ' (1) 

where /bk is the reduced bulk free energy per volume, while / S)0 + / S; b is the surface excess free 
energy that results as the separation L of the two confining plates is increased to infinity. According 
to the theory of finite-size scaling, the residual free energy should take the scaling form 

fl%L) « L-^- 1 ) QadL/ioo) (2) 

on sufficiently large length scales, where ^ is the bulk correlation length. Here Q a .b is a universal 
scaling function, which depends on the bulk universality class, the geometry of the system, and gross 
properties of the confining planes, such as large-scale boundary conditions (BC) associated with RG 
fixed points of the corresponding boundary field theory, but is independent of microscopic details. Its 
value at bulk criticality, A 0) & = O a ,b(0), is a universal number, the so-called Casimir amplitude. 

The so far most comprehensive and detailed theoretical investigation of critical Casimir forces is 
that of Krech and Dietrich (KD) [8, 9] who considered n-component cf) 4 models in a oo d_1 x L slab 
geometry whose Hamiltonian is given by 



H[4>] = [ dz [ d^r 
Jo if 1 - 1 



2^ +2^ + 4! 



(3) 



together with BC. KD considered periodic (per) and antiperiodic (ap) BC, as well as (a,b) = (D,D), 
(D,sp), and (sp, sp). The notations D and sp indicate Dirichlet and special BC on the boundary 
planes Si at z = and *B 2 at z = L. The Dirichlet BC corresponds to the infrared-stable fixed 
point describing the surface critical behavior at the ordinary transition of semi-infinite systems. By 
special, the large-scale BC on a plate is meant that applies when the surface interactions on it are 
critically enhanced and no symmetry breaking surface terms are present there; they pertain to the fixed 
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point associated with the so-called special transition [22]. All combinations (a, b) of such BC with 
a, b = D, sp can be implemented by adding surface terms J2j=i I<s 2 /2 to the Hamiltonian 
with tj — oo or Cj = c sp (critical enhancement [25]), depending on whether a, b = D, sp. 

Performing two-loop calculations of the excess free energies, Casimir forces, and their scaling 
functions for the BC p = per, ap, (D, D), (D, sp), and (sp, sp), and restricting themselves to tem- 
peratures T > T c ,oo, KD obtained expansions to first order in e of the Casimir amplitudes A p and the 
scaling functions 9 p (y). 

The restriction T > T c oo had not only technical reasons. As T is lowered beneath T c oo , a 
crossover to the critical behavior of an effectively d — 1 dimensional system is expected to occur at 
a temperature T c l < T c>0o . The e expansion is unable to deal with the d — 1 dimensional infrared 
singularities at T Cj l since the appropriate small dimensional parameter would be €5 = 5 — d. KD 
were aware of this problem. They verified that their 0(e) results for the boundary conditions p = 
ap, (D, D), and (D, sp) were consistent with the asymptotic behavior of the scaling functions Q p (y) 
for y — > 0± one can infer from the requirement that fi be analytic at T c oo when L < 00. They also 
noted that one could not expect their results for per and sp . S p to exhibit such a small-?/ behavior 
because the corresponding free propagators with periodic and Neumann boundary conditions could 
not be analytically continued to negative values of r ~ (T — T c ,oa) /T c>00 , unlike those for the other 
boundary conditions. In order to assess the predictive power of their 0(e) results for P , KD (and 
others [12,26]) nevertheless extrapolated these functions for all five types of boundary conditions to 
d = 3 by setting e = 1. 

To see that the e expansion breaks down at T = T c-OQ for p — per and (sp, sp), consider the 
three-loop graph OOO °f /i> where the lines represent the free propagator 

between the points x = (r, z) and x' = (r', z'). Here (z\m) and are the orthonormal eigen- 
functions and eigenvalues for the BC p, respectively. For example, for p = per one has (z\m) = 
£-i/2 e ifc m 2 with k m = 2-kiji/L and m £ Z, whereas (z\m) = y/(2 — 5 m ,a)/L cos(fc m z) for 
p = (sp, sp), with k m = nm/L, m = 0, 1, . . . , oo. 

The central part of this three-loop graph is the subgraph ^y, which behaves in the bulk case 
L = 00 (for zero external momenta) as f -E / 2 , and hence is infrared singular at T c oo [27]. However, 
each of the two other bubbles contributes a factor G^°°^(a;; x) ~ f 1_e / 2 . Hence, the limit f — ► 
of the overall contribution vanishes [27]. When L < 00, each of the two tadpoles G^\x;x) for 
p = per, (sp, sp) can be decomposed into a contribution (PoGp Pq)(x; x) ~ f^ 1_e ^ 2 from the 

fe m = mode and a remainder (QoGpQo)(x; x), where Pq = 1 — Qo = |0)(0|. At f = we are 
left with the contributions from the remainder, namely 

(QaG£lQo)(x;x)k= = ^"ff (5) 

and 



(QoGW p Q )(x; x)|, =0 = ^i-J^ [ 2 C (2 - e) + c(2 - e, J) + c(2 - e, ^) 

where C(s, a) = YljLoU + ^ i s tri e Hurwitz zeta function. 

The zero-mode contribution [(PqG^ Po){x; x')] 2 of ^3 must be integrated over the parallel sep- 
arationr — r'. The result behaves ~ T~( 1 + t )/ 2 . Combined with the (QoGpQo) parts of the tadpoles, 



(6) 
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it produces a contribution that diverges as f — > 0. Thus the loop expansion is ill-defined at f = for 
these two BC p = per and (sp, sp). Evidently, this breakdown should occur more generally whenever 
ak m — mode is present, which is not the case for p — ap, (D, D), and (D, sp). 

The origin of this breakdown is a deficiency of Landau theory: Whenever the free propagator 
involves a k m = mode, it predicts a transition for the bulk and the film of finite thickness L at the 
same critical value f — 0. The remedy is a reorganization of the perturbation series. The k m = mode 
must be split off and treated in the background of the k ^ modes. As a result of its coupling to the 
latter, the zero mode becomes massive at T c oo for L < oo. 

To formulate such an expansion, we decompose <f> as <p{x) = <p(r) 
contribution p(r) and its orthogonal component tp with f„ L dz if)(r, z) - 
us an effective d — 1 dimensional field theory with the Hamiltonian 

7Us[<p] =-mTi >e - H ^]. 

Let us introduce the free-energy part and the average (. . )^ by 



f ip(x) into its k = mode 
0. Integrating out tp defines 

(7) 



F^ = - lnTr^,e 



Ti>(. 



Then we have 



with 



H eS [tp] = F^ + H[<p] - \n{e- HM ^)ip 



and 



H[ V ] 



dz 



G 



L 



tL 



uL 
4! 



<P 



(8) 



(9) 



(10) 



(11) 



The last term in Eq. (|9j gives loop corrections Y^bLi ^iff 
obtains 



For the one-loop contribution, one 



Hi 



= -Trl n [l + (n/6)G^(<W 



2<y9 a (^j)] 



I I 



(12) 



where the dashed lines represent free V'-propagators G^' — QoGpQo, while the gray bars indicate 
•p legs. 



Added to %[<p], the first graph in Eq. ill 2b produces the shift 



t { l) = t + 5fi L) with 



so that the free (^-propagator G^ 1 acquires an i-dependent mass at f = 0. The second graph in 
Eq. d!2t yields a nonlocal p 2 p 2 interaction; the suppressed ones correspond to similar nonlocal inter- 
actions involving more than two p 2 operators. The two-loop term Hgg[^] involves two contributions 
~ tpip; a local one giving an 0(u 2 ) correction to the shift 8fp, and a nonlocal one whose interaction 
potential is proportional to J Q dz J Q dz' [G^(x,x')) 3 . 

Upon substituting the results and ijfji into Eq. i\3l . the shifts Sfpev and 67 S p,sp can be computed 
in a straightforward manner. At f — 0, one obtains 



P 



u n 
2 _ 



I -r {QoG v )Qo){x;xh (13) 



"'per z °'sp,sp — u 



2r(l-e/2)C(2-e) 



6 



27r 2- e /2 L 2- e 



(14) 
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We can now set up a Feynman graph expansion, utilizing G v {p) = \L(p 2 + Tp^)] 1 as free 
propagator in the d — 1 dimensional momentum space and employing dimensional regularization. 
Since G v {p) remains massive at r = when L < oo, the Feynman graphs are infrared finite as 
long as f > 0. Owing to our reorganization of perturbation theory, the breakdown encountered in the 
conventional expansion in terms of G p is avoided. Note also that the correct bulk expressions are 
recovered as L — ► oo since the contributions from the zero mode vanish. For example, the term in 
Eq. (|9j yields the bulk free energy when L — > oo. Indeed, setting T = T C)00 and denoting the analog 
of f p es for Fj, as f$? p we find 



L 



d-l 



i?res 

'i>;p\T r 



,(0) 



(e) + ^a«( e )to« + 0(* a £*) 



(15) 



with 



a (o) (e) = 2 4 - £ a (°) (r) - r ^ 2 ~ 6 ^ 

"perW " "sp,spl c /> — _.2— e/2 

2/ 



,(1) 



,(1) 



(e) = 
,(*) = 



Z 7T 



p,sp 
£-4 t^2 



C(4-e) : 



T 2 (l - e/2)C (2 
2C 2 (e-l) + C(2e 



(16) 



4 3 - £ 2cos 2 (7re/2)r 2 [(3-e)/2] 

KD's two-loop results for f p \r a with p = per, (sp, sp) follow from Eqs. (1151 and (1161 . as they 
should because the zero-mode contributions to f™ s of order u° and u vanish at T c oo . However, the 
^-dependent part of H c s[(p] gives additional contributions. The leading ones correspond to one- and 
two-loop terms of a <ys 4 theory in d — 1 dimensions with mass coefficient fp and coupling constant 
u/L, as the rescaling L x l 2 Lp — > in Eq. (II Q shows. At T^qo, they yield L-dependent contributions 
to /™ s proportional to (<w£ £) )( 3 - e )/ 2 - L^^uL^ 3 - 6 )/ 2 and (u/L)((5f f i L) ) 1 - £ - L 1_d {uL e ) 2 ~\ 
respectively. Including only the first one, we arrive at 

L d - 



r rvcs _ /res 1 
'Up ~ Jtp;p\T c 



A P {e) 



4! 



-uV 



(3-£)/2 



with 



A pcr (e) = 2< 2 - 



)(M/2j c e i - _ 

jri sp,spv c / 



r[(e-3)/2] 

2 4~e 



2r(l-e/2)C(2-e) 



7r (6-e)/2 



(3-£)/2 



(17) 



(18) 



To combine these results with RG-improved perturbation theory, we follow KD. We utilize the 
reparametrizations (u = 2 d n d / 2 Z u fj, e u, f = [x 2 Z t t, . ..) of bulk and surface quantities of the cor- 
responding semi-infinite theories [22], and fix the additional additive (bulk and surface) counterterms 
Jl requires such that /i )re n> its renormalized counterpart, vanishes at r = 1 together with its 1st 
and 2nd r-derivatives [28]. To obtain the critical Casimir amplitudes A p , we must express /™ s in 
terms of renormalized variables, set [ih = 1 and t = 0, and evaluate it at the fixed-point value 
u* = 3e/ (n + 8) + 0(e 2 ). Upon expanding in powers of e, we find (7 = Euler-Mascheroni constant) 



90 



180 



ln7r 



2C'(4) + 5n+2 



2n + 8 



ir 2 f 1 
9VG V 



3/2 



p 3/2 



0(e 2 ) (19) 



and 



A 



sp,sp 



7T 



n 2 e 



1440 

_2 



2880 



1 - 7 - ln(47r) + - 



2C'(4) 



2n + 8 C(4) 



7T" 



72 V6 



3/2 



,3/2 



0(e 2 ) 



(20) 
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The 0(e 3 / 2 ) terms result from the 0(iS z ^f 2 ) contributions in Eq. d!7l >. Obviously, the latter also 
implies contributions of the form e fc + 3 / 2 ln fc e with fceN. Furthermore, the terms (u/L)(5fp i ' ) ) 1 ~ £ ~ 
u 2 ~ t mentioned above yield contributions of the form e k+2 ln fe e. 

Can the appearance of the e 3//2 and unconventional higher-order terms be checked by alternative 
means? This is indeed possible: The limiting value lirrin^oo A per /n can be obtained from the ex- 
act solution of the mean-spherical model (expressed in terms of the function Yq of [29], it becomes 
Yo(0, 0)). Solving the corresponding self-consistent equations iteratively with e > reproduces the 
n — ► oo limit of all terms on the right-hand side of Eq. d!9l > and shows the existence of higher-order 
contributions of the mentioned form [30]. 

In Table|I|we give the values of A p Eqs. dl9l and J20i predict for n = 1, 2, 3, oo upon setting e = 1. 

TABLE I - Estimated values of A p (d,n) /n for d — 3. Following [9], we have included the d — 3 values of 
A sPiSp , although they are probably of little physical interest because the special surface transition is not expected 
to occur at d = 3 in the 0(n > 2) case, unless surface couplings of infinite strengths are allowed. 



n 


1 


2 


3 


oo 


Ap CT (3, n)/n 


-0.1967" 


-0.2147" 


-0.2311" 


-0.4668" 




-0.1105 6 


-0.1014 6 


-0.0939 6 


-0.0192 6 




-0.1526 c 






-0.1531 d 


A sp , sp (3,n)/n 


-0.0224 a 


-0.0252 a 


-0.0278" 


-0.0619" 




-0.0117 6 


-O.Olll 6 


-0.0106 6 


-0.0059 6 



" Values obtained by setting e = 1 in Eqs. i 191 and I20i . 

b 0(t) results [9], evaluated at e = 1. 
c MC results according to [31]. 

d Exact value -2C(3)/(5tt) according to [32] and [29]. 



For comparison, the corresponding O(e) estimates are also listed, along with a MC estimate [31] and 
an exact n = oo result [29,32]. A known problem of the 0(e) results for A pcr is the seemingly 
incorrect n-dependence of the predicted d = 3 values, whose deviations from the exact n = oo value 
increase monotonically as n grows (see Fig. 12.8 of [12] and [16]), although the MC estimate for 
n = 1 is very close to the exact n = oo value. The e 3 / 2 term modifies the n dependence, yielding an 
estimate for — A pcr /n that increases with n. 

In summary, we have shown the following: (i) the e expansions of quantities such as Casimir ampli- 
tudes are ill-defined at T c oo when the BC gives a zero mode in Landau theory, (ii) The reformulation 
of field theory presented here yields well-defined small- e expansions for temperatures T > T Cj00 . In 
the cases of A pcr and A sp sp , these expansions involve fractional powers and logarithms of e. Clearly, 
more work is necessary to explore the potential of such expansions for reliable extrapolations to d = 3. 

In typical experimental situations one expects to have either Robin BC d n <p = Cj<p (which in 
the long-scale limit normally should map on (D, D) BC) or symmetry-breaking (+, ±) BC (classi- 
cal liquids). However, experimental situations corresponding to near-critical enhancement of surface 
interactions on both plates are conceivable. In that case a crossover from an initial behavior character- 
istic of (sp, sp) BC should occur. Clearly, proper treatments of this crossover must take into account 
the findings described above. Needless to say, that in MC simulations periodic BC are the preferred 
choice, and that simulations dealing with (sp, sp) BC were performed as well [33]. 

* * * 
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